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Abstract 

An unbroken Z 3 symmetry remains when a local SU(2)x symmetry is broken spontaneously by 
a quadruplet. The gauge boson Xh(X/j.) carries the dark charge and is the candidate of dark matter 
(DM). Due to the mixture of the scalar boson <f> r of the quadruplet and the standard model (SM) 
Higgs boson, the DM can annihilate into SM particles through the Higgs portal. To investigate the 
implications of the vector DM in the model, we study the relic density of DM, the direct detection 
of the DM-nucleon scattering and the excess of the gamma-ray spectrum from the Galactic Center, 
which is supported by the data from the Fermi Gamma-Ray Space Telescope. We find that with 
the DM mass of around 70 GeV in our model, the data for the excess of the gamma-ray could be 
fitted well. 
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One of unsolved problems in astrophysics is the existence of dark matter (DM), where 
the plausible candidates in particle physics are the weakly interacting massive particles 
(WIMPs). The Planck best-fit for the DM density, which combines the data of the WMAP 
polarization at low multipoles, liigh-f' experiments and baryon acoustic oscillations (BAO), 
is given by [l| 

n DM h 2 = 0.1187 ±0.0017. (1) 


Besides the evidence from astronomical observations, now there are direct and indirec 


ways to detect DM. According to the recent measurements by XENONIOO [2f] and LUX [ 3 ] 
Collaborations, which are designed for directly detecting DM, since no clear signal is found, 
the cross section for the elastic scattering of DM off nucleons has been strictly limited. 
Additionally, although the potential DM signals are indicated by the indirect detections, 
such as the excess of the positron fraction observed by PAMELA d| and Fermi-LAT j(j| 
experiments, and the excess of the positron+electron flux observed by PAMELA |7], Fermi- 
LAT {sl-ATIC [9 . 
pulsars [12, 13]. 


and HESS 


10 


Ilf . they may also be solved by astronomical effects, e.g. 


photon energy of around 2 GeV, has been pointed out by the analyses in Refs. 


at 

the 

14 

-21]. 


Furthermore, using the data from the observation of the Fermi Gamma-Ray Space Tele¬ 


scope 


22 


23], a more significant signal of the gamma-ray from the region around the Galac¬ 


tic Center is also found 


24 


,28| . Subsequently, it has been found that the excess matches 


well with the gamma-ray spectrum from the DM annihilation, where the requested thermally 
averaged cross section (<7f re i) at the order of 10 -26 cm 3 /s is the same as that of the thermal 
relic density. Moreover, it has been pointed out that the effects through the Higgs portal 
could naturally explain the excess of the gamma-ray spectrum 24], l37j. Based on these 


results, in this paper we propose a stable vector DM model in which a discrete symmetry 
stabilizing the DM is obtained naturally and the DMs annihilate into SM fermions through 
the Higgs portal. Other mechanisms to explain the excess could be referred to the references 
in the literature, such as that DM annihilates directly into SM particles and/or DM first 
annihilates into hidden scalar (gauge) bosons, and then decays to SM particles via the Higgs 
(Z')-portal 29 56]. 

From the view point of model buildings, to protect DM from its decay, an unbroken 
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symmetry in the theory is necessary. However, a discrete symmetry usually is put in by 
hand. In order to get a stable DM naturally, we study the model in which the unbroken 
symmetry originates from a spontaneously broken gauge symmetry. To realize the concept, 
particularly we are interested in the extension of the SM with a new SU(2)x gauge symmetry 
where the subscript X is regarded as a dark charge. The interesting properties of a local 
SU{2)x group are: (1) comparing with the local 17(1) case in which the U(l) charge has to 


satisfy some artificial tuning 


57], an unbroken discrete symmetry can be naturally preserved 


after the spontaneous breaking of the SU(2)x gauge symmetry; (2) the massive gauge bosons 
from SU(2) X could be the DM candidates. The various applications of the hidden SU(2 ) 


rave been studied in the literature, sue 


58|, a custodial symmetry in Refs. 


32 


l as a remaining Z 2 symmetry with 


59 ] and an unbroken 7/(1) of SU(2 ) 


gauge symmetry 
a quintet in Ref. 
in Refs. 60, 61]]. 

Since the model with the custodial symmetry discussed in Ref. |59j is similar to our 
proposal, it is worthy to show the difference between them. It has been noticed that without 
introducing any new fermions or higher multiple states in the hidden SU(2)x gauge sector, 
a new fundamental representation of SU (2 W co uld lead to three degenerate DM candidates 
by utilizing the SO (3) custodial symmetry |59j. Due to the custodial symmetry, the three 
DM candidates are stable particles. However, the symmetry could be broken easily when 
SU(2)x fermions and/or higher representation scalar fields are included. Although the 
inclusion of the new fermionic and/or higher multiple staff is not necessary, if one connects 
the origin of neutrino masses with the dark sector, the inclusion of the new staff becomes a 
relevant issue. In order to get over the possible unstable effects when more phenomenological 
problems in particle physics are involved, we propose to use a discrete symmetry to stabilize 
DM, where the discrete symmetry is not broken by higher multiplet fields or fermions under 
SU(2)x- Additionally, the processes for explaining the gamma-ray excess in our model are 
different from those dictated by the custodial symmetry 32, [59]. We will see the differences 
in the analysis below. Moreover, we find that an Z 3 discrete symmetry indeed remains when 
SU(2)x is broken by a scalar quadruplet. Based on the introduced quadruplet, we summarize 
the characteristics of our model as follows: (a) the unbroken Z 3 symmetry is the remnant 
of SU(2)xi (b) two gauge bosons Xu and Xu carry the Z 3 charge and are the candidates of 
DM, (c) besides the SM Higgs (</>), only one new scalar boson (<f> r ) is introduced, and (d) 
due to the mixture of 0 r and (f), the DM annihilation is through the Higgs portal. 
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In the following, we briefly introduce the model and discuss the relevant interactions with 
the candidates of DM. To study the minimal extension of the SM that includes the staff of 
DM, besides the SM particles and their dictated gauge symmetry, we consider a new local 
SU(2) X gauge symmetry and add one quadruplet of SU(2) X to the model. The introduced 
quadruplet is not only responsible for the breaking of the new gauge symmetry, but also 
plays an important role on the communication between dark and visible sectors. Thus, the 
Lagrangian in SU(2) X x SU(2)l x f/(l)y is written as 

C = C S M + (D^ 4 - V(H, $ 4 ) - \x%X a ^ (2) 

with 


V(H, $ 4 ) = n 2 H*H + \(H*H) 2 + + A $ ($Id > 4 ) 2 + A ,< f>4<3>4if t if, (3) 


where C S m is the Lagrangian of the SM, H T = (G + , (v + 0 + iG°)/\/2) is the SM Higgs 
doublet, $4 = (03/2,01/2, —0-1/2, 0 - 3 / 2 )/v / 2 is the quadruplet of SU(2) X , the index i of 0 * 
stands for the eigenvalue of the third generator of SU{2) x , 0-i = 0*, the covariant derivative 
of <h 4 is Dn = dp, + ig x T a X “ with the representations of T a in the quadraplet, given by 
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and T 3 = diag(3/2,1/2, —1/2, — 3/2), and the field strength tensor of SU(2) X is read by 

x; u = d,x a u - d„x; - gx (x, x x v y. 

To break SU(2) X but preserve a discrete symmetry, the non-vanishing vacuum expecta¬ 
tion value (VEV) and the associated fields fluctuated around the VEV are set to be 


(0±3/2) — , 0±3/2 — + 0r ± l£) ■ 

When we regard the quadruplet as the fluctuations from the vacuum $0 
$ 4 can be parametrized by using the form 


(5) 

(n 4 , 0, 0, w 4 )/2, 


$ 4 = e « ra ofi(x)/v*Q 4 ; ( 6 ) 

$4 = -^= (0 r ,O,O,0 r ) 
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with cf) r = (i>4 + (f> r )/\/ 2. In terms of scalar fields a a (x), the components of <h 4 could be 
expressed as </>i /2 = \/3(— a 2 (x ) + ia 1 (x))/ 2-^/2, (j> - 1/2 = and £ = 3/2a 3 (a:), where we 
have taken the leading terms in the held expansions. Eq. (j6]) indeed is nothing but a local 
gauge transformation. Therefore, </>±i/2 and £ could be rotated away from the kinetic term 
of $4 and the scalar potential; and they are the unphysical Nambu-Goldstone (NG) bosons 
of the local SU(2)x symmetry breaking. Consequently, we can just employ S> 4 for exploring 
the mass spectra of new particles. 

With the breaking pattern in Eq. ([5]), one can find that an Z 3 symmetry U 3 = e* T34,r / 3 = 
diag(l, e* 27r / 3 , e~ 2l7r / 3 ,1) is preserved by the ground state <h 0 . Under the Z 3 transformation, 
the scalar fields of the quadruplet are transformed as 


0±3/2 -> 0±3/2 , 

0±i/2 —► e ±i2 " /3 0±i/2 • (7) 


That is, 0±3/2 are Z 3 blind while (ft± 1 carry the charges of Z 3 . To understand the transfor¬ 
mations of gauge fields, one can use 

T a X'“ = U 3 T b X*ul. (8) 

In terms of physical states of gauge fields, one can write 

T“x; = 4(T + X„ + T~x r ) + T 3 Xl (9) 

with T± = T 1 ± iT 2 and Xm(Xm) = -F iX 2 )/\/2 where is regarded as the antiparticle 
of Xfi- Using the identity U^T^l = exp(±f47r/3)T ± , the transformations of X^iX/i) and X 3 
under Z 3 are given by 


^/x ' ^/x ? 

X/Xw) —* e±i4n ^ 3 Xn(Xn) • (10) 


We see that X/j,(Xh) carries the Z 3 charge and X 3 is the Z 3 blind. Due to the unbroken Z 3 , 
the particles with the charges of Z 3 are the candidates of DM. Since <fi±i /2 are the unphysical 
NG bosons, the DM candidates in our model are the vector gauge bosons Xu and 

To study the spectra of SU(2)x, we have to determine the nonvanishing VEVs of H and 
<h 4 . Using Eqs. (EJ) and (jSj), we get 


, s. v 2 n 2 Xv 4 n%v\ X^v\ \'v 2 vl 
V(v,v*) = ^— + — + ^—+ 


( 11 ) 


5 







With minimal conditions dV(v,V 4 )/dv = dV(v,V 4 )/dv 4 = 0, we have 


/i 2 + Xv 2 + 
+ A$W 2 + 


X'v'i 

2 

X'v 2 


o, 

0, 


( 12 ) 


respectively. In terms of the parameters in the scalar potential, the VEVs could be written 
as 


2A'/i| — 4A$/i 2 
4AA$ - A' 2 
2X'/li 2 - 4A/i| 
4AA$ — A' 2 ' 


(13) 


As known that the masses of gauge bosons arise from the kinetic term of $ 4 , accordingly 
the masses of Xn{Xn) and A 2 can be directly found by 


%g 2 x [(T~T + + T + T-)x»r + (T 3 ) 2 X 3 X 3 ^] $0 


!)Wa 


2(t(t + 1) - ti)x,r + tlxlx^ 


(14) 


where t(t + 1) and t 3 are the eigenvalues of T 2 = T a T a and T 3 , respectively. With t — t, :i — 
3/2, the masses of gauge bosons are obtained as 


V3 


"A — ~^9xV4 , rn X 3 — -gxV4 


(15) 


Although there are four scalar fields in the quadruplet, three of them become the longi¬ 
tudinal polarizations of gauge bosons (y M , W 3 ). Therefore, combining with the Higgs 
doublet in the SM, the remaining physical scalar bosons in the model are 0 and (j) r . In terms 
of the scalar potential in Eq. ([3]), the mass matrix for 0 and <p r is expressed by 


M 2 


( m 2 X'vv 4 
X'vv 4 m\ r 


(16) 


with = y/2Xv and = yj 2A$u 4 . Due to the X 1 effect, the SM Higgs 0 and 0 r will 
mix and are not physical eigenstates. The mixing angle connected with the mass eigenstates 
could be parametrized by 



cos 9 sin 6 
— sin 9 cos 9 



( 17 ) 
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where h denotes the SM-like Higgs, H° is the second scalar boson and tan 26 = 2\'vv±/{m\ — 
m ^). According to Eq. (jl 6 ]l . the mass squares of physical scalars are found by 

m l,2 = \{™l + m lr ± + AX ,2 V 2 VI) . (18) 

We note that the mass of h could be mi or m 2 and the mass assignment depends on the 
chosen scheme of the parameters. To solve the problem of the gamma-ray excess, we will 
focus on the case of > m H 0 . 

Next, we derive the couplings of <fi and (j) r and the interactions with new gauge bosons. We 
first discuss the gauge interactions of <j> r . From Eq. (J2]), we see that the gauge interactions 
of the quadruplet only occur in the kinetic term of $ 4 . Using 3>4 defined in Eq. (J 6 J) and the 
covariant derivative of $ 4 , the gauge interactions are expressed as 

Ig = dfi\ (igT a X atl ) <F 4 + h.c ., (19) 

I go = {tgT a X;$ 4 ) ] (igT b X b ^ 4 ) . (20) 


By adopting the expression of Eq. (JUJ), one can easily find that the gauge interactions of 
Eq. (fl9|i vanish. By using the result 


T a X“$ 4 


3/2 xl \ 


-3/2 Xl) 


V 4 + 
2 


(j) r 


? 


( 21 ) 


Eq. flUD can be straightforwardly written as 

3 yd) 

Igg = VSgxm^rXnX 1 * + -^-axm^rX^X 3 ^ 

+ 5 (^) <Hx„r + \ (^) . ( 22 ) 

where the masses of gauge bosons defined in Eq. (1T5l) have been applied. We second discuss 
the couplings of <f) r to the SM Higgs 0 where the vertices could be obtained from the scalar 
potential of Eq. ()3|) . Since the derivations are straightforward, we summarize the vertices 
of (f) r and (j) in Table Q] We note that although the interactions in Eq. (122|) and Table U are 
shown in terms of <f> r and 0 , the expressions with h and H° mass eigenstates could be easily 
obtained when Eq. (fT7ji is applied. 
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TABLE I: Couplings of the scalar boson <fi r to SM Higgs <fi. 


<\>A 2 


CO fc. 

-e- 



AA 4 

X'v 

3!A$U4 

A' 

3!A$ 


The relevant free parameters in the model are /ijU, A($), A' and the gauge coupling g x . 
Using the masses of </> and </y and the VEVs of H and $ 4 , the six parameters could be 
replaced by (g x , v, tq, m^, m^, A'). When these values of parameters are fixed, the masses 
of h and H° and the mixing angle 6 are determined. According to the results measured by 


ATLAS 


62| and CMS 63], the Higgs mass now is known to be rrih = 125 GeV. Therefore, 


it is better to use the physical masses m^H 0 and mixing angle 9 instead of and A'. 

Additionally, the VEV of v ~ 246 GeV is determined from the Fermi constant Gp and tq 
can be replaced by m x . Hence, the involving unknown parameters in the model are g x , m x , 
mp 0 and 6. 


To constrain t 
the relic density 


re free parameters, two observables have to be taken into account: one is 
l| and another one is the DM-nucleon scattering cross section 2, 3]. The 


number density of DM is dictated by the well-known Boltzmann equation, expressed by 


(b + ally, — -(/r// re i) ( n 2 - nl q ) 


(23) 


where H is the Hubble parameter, n — n x + n x , and n eq is the equilibrium density, defined 
by 


n X,eq — n x,eq — 9x 


m\T 

~ 2^ K2 


m 


(24) 


with g x the internal degrees of freedom of DM, T the temperature and AA the modified 


Bessel function of the second kind 


64j . For the vector DM, we take g x = 3. The thermally 


averaged annihilation cross section is given by 

1 < 


(crV Te 1 ) = 


8 Tm\Kl(m x /T) J 4ml 


dsy/s(s - Am x )Ki(y/s/T)a(xX all). 


(25) 


In the model, the DM annihilating into the SM particles is through the Higgs portal, where 
the associated Feynman diagrams are presented in Fig. [T] We note that in contrast to 


Ref. 


32], the DM semi-annihilation processes such as XX x{H°, h ) are absent in our 


model. To study the DM abundance after the freeze-out, usually it is more convenient to 



















FIG. 1: Processes of the DM annihilation. Diagram II includes t- and u-channel. 


consider the ratio of the number density to entropy density, defined by Y = n/s, where 
s = (27r 2 /45)g*T 3 and g*(T ) is the effective number of degrees of freedom contributing to 
the entropy density. With H = — T/T, s + 3 Hs = 0 and x = m x /T, Eq. Q23J) leads to 


dY_ 

dx 


47 r m x Mp 
\/90 x 2 


V^XT){av„ i) (Y 2 



(26) 


where H 2 = 87r 3 Gg*T 4 /90 and Mp = 1/(87 tG) have been used. If we set Y 00 to be the 
present value after the freeze-out, the current relic density of DM is given by 


q _ ^xi x SqY 00 

x ~ 3 H 2 M 2 ’ 


(27) 


where H$ and So are the present Hubble constant and entropy density, respectively. For 
numerical calculations, we employ micrOMEGAs 4.1.5 65| to solve the Boltzmann equation 

and get the present relic density of DM defined in Eq. (12711 . 

Although the direct detection of DM via the DM-nucleon scattering has not been observed 
yet, the sensitivity of the current experiment could give a strict constraint on the free 
parameters. In the model, the sketch of a vector DM scattering off a nucleon is shown 
in Fig. [2j By neglecting the small momentum transfer, the scattering amplitude of the 
X /i (x M )-nucleon is written as 


M = 


V3g x m x m 2 h - m 


2 

H° 


rn l m H o 


sin 9 cos 9(N\m q qq\N) 


(28) 


By assuming that the effective couplings of DM to the proton and neutron are the same, 


we parametrize the nucleon transition matrix e 
where the range of is [1.1, 3.2] x 10 


-3 


66 


ement to be ( N\m q qq\N) = f N / f) 1 ^ 2 , 
67]. As a result, the scattering cross section of 
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FIG. 2: The sketch of the dark matter scattering off a nucleon. 


the DM-nucleon is formulated by 

3 g 2 x f 2 


°XN ~ °~X(X)N^ X (X)N 


N 


sin 6 cos 6)~ 


m N 


m H o ~ m~ h 


(29) 


47 T ' ' \m x + rriNj \ m h mHO 

Before discussing the numerical analysis, we set up the possible schemes for the values 
of m x and o. Since XX W + W~, ZZ are the dominant channels in the case of m#o > 


m x > mw and in disfavor with the gamma-ray spectrum 


28] . we assume x is lighter than 


W and Z. To explain the excess of the gamma-ray spectrum, it has been pointed out that 


the preferred channels via t 


scalar and yy y bb [28|, 


32 


re Higgs portal are XX —> SS —> bbbb with S being the possible 


37], where the former produces the on-shell S and subsequently 


S decays into SM particles while the latter utilizes the resonant enhancement of ms ~ 2 m x . 
As a result, we focus on the following two schemes: 

(a) m x = 70,60 GeV and m^o < m x , where the DM annihilation channel is XX H°H° 
with H° being the on-shell scalar boson; and afterwards H° decays through H° —> bb [ 32 , 37]. 
k} HZx = 50,40 GeV and o > m x , where the DM annihilation channel is XX 


bb 


m 


32]. We will see that the channel becomes significant when the condition of o 


2 m x is satished. 

Although the fermions in the final states could be other lighter leptons and quarks, since 
the coupling of the scalar to the fermion depends on the mass of the fermion, we only focus 
on the 6-quark pairs in the final states. 

In scheme (a), as the main DM annihilating processes are from Figs. [3 and [TJI and 
the produced H° pairs are on-shell, the results are insensitive to the mixing angle 6. To 
understand the constraint of the observed Ddm^ 2 , we present Q x h 2 as a function of gx i 


m 


Fig. [3](a). From the results, we see that for matching the observed relic density of DM, 
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the value of the gauge coupling g x should be around 0.23(0.21) for m x = 70(60) GeV and 
run o = 69(59) GeV. We note that for explaining the excess of the gamma-ray via the DM 
annihilation, we adopt m x ~ m H o in scheme (a). We will clarify this point later. In scheme 
(b), Fig. [UII becomes dominant. Since h and H° both contribute to the DM annihilation, 
besides the gauge coupling g x and m H o, the results are also sensitive to the mixing angle 
6. Since there are three free parameters involved in this scheme, in Fig. |3](b) we show 
the correlation between sin 6 and m#o when g x = 1 is taken and the observed Ddm^ 2 is 
simultaneously satisfied. 



gx m H o [GeV] 


FIG. 3: (a) Relic density of x(x) as a function of g x in scheme (a), where the band indicates the 

observed value of Ddm^ 2 - (b) Correlation between sin@ and m H o in scheme (b) when g x = 1 is 
taken and the observed relic density of DM is satisfied. 


With the proposed schemes (a) and (b), we can further discuss the constraints from the 
measurements of DM direct detections. Since the vector DM candidates are not self-charge- 
conjugation particles, the DM density is composed of Xu and Xm be. pdm = p x + p x . Thus, 
the elastic scattering cross section of DM off a nucleon is proportional to P x ct x n + Px a xN — 
PdmO'xn- Consequently, for comparing with the DM-nucleon scattering cross section mea¬ 
sured by the direct detection experiments, one can just use a xN which is formulated in 
Eq. (I29lh For scheme (a), unlike f\h 2 , a xX is sin 0 dependent. We plot the elastic cross 
section as a function of sin 0 in Fig. |4j where we have taken (m x ,m#o) = (70,69) GeV for 
the left panel and (60, 59) GeV for the right panel. In order to fit the measurement of Ddm h 2 
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simultaneously, we use g x = 0.23(0.21) for the former (latter). For comparisons, we also 


show the 90%-CL upper limits by XENONIOO |2j and LUX [3] Collaborations on the plots. 
From the results, we clearly see that to satisfy the DM direct detection experiments, we need 
sin 0 <0.1. For scheme (b), we present a x x as a function of m H o in Fig. 0 with g x = 1 and 
m x = 50 (40) GeV for the left (right) panel. In order to fit the data of U DM h 2 together, in 
the figure we have applied the results shown in Fig. [31(b). By the plots, we find that current 
DM direct detection experiments further limit the mass relation to be m#o ~ 2 m x 


"x- 



FIG. 4: DM-nucleon scattering cross section in scheme (a), where the constraint of observed 
Ddm^ 2 has been considered. For comparisons, the measurements of XENONIOO [2] and LUX 
for 90%-CL upper limits are shown in the plots. 


After analyzing the constraints of the DM relic density and direct detection, we now 
study the gamma-ray which is originated from the DM annihilation. It is known that the 
flux of the gamma-ray from the DM annihilation is expressed by 

d$(E 1 ,'ip) (av ie \) dNry r 


dE^dVl 


8vrm^ dE x J los 


P (r)dl(ip ), 


(30) 


where dN x /dE x is the gamma-ray spectrum produced per annihilation, ip is the observation 
angle between the line-of-sight and the galactic center, p(r) is density of DM, and the 
integration of the density squared is carried out over the line-of-sight. The general DM halo 
profile could be parametrized by 


^ (r®yt (1 + (r Q /r s ) c 

pM = p©Lt) (ttfuf 


q \ 03-7)/“ 


(31) 
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FIG. 5: DM-nucleon scattering cross section as a function of m^o in scheme (b), where the con¬ 
straint of observed Udm h 2 has been considered. The measurements of XENONIOO [2] and LUX 
for 90%-CL upper limits at the corresponding m x are also shown in the plots. 


where r s = 20 kpc is the scale radius, p & = 0.3 GeV/cm 3 is the local dark matter density at 
r Q = 8.5 kpc and r is the distance from the center of the galaxy. Note that (a, (3, 7) = (1, 3,1) 
corresponds to the Navarro-Frenk-White (NFW) profile. In our numerical estimations, we 
set a = 1 and j3 = 3, but 7 to be a free parameter. Since p[r) is proportional to r 


-7 


we 


see that the change of the parameter 7 can only shift the entire gamma-ray spectrum but 
not the shape of gamma-ray flux. For executing the numerical calculations of Eq. ([30]) . we 
implement our model to raicrOMEGAs 4.1.5 65] and use the program code to estimate 

the gamma-ray spectrum. 

In the model, the processes to produce the gamma-ray by the DM annihilation are similar 
to those for the relic density, except that the gamma-ray is emitted in the final states. In 
scheme (a), we present the flux of the gamma-ray as a function of the photon energy E 1 in 
Fig. [6j]a) , where the solid line denotes (m x ,m H 0) = (70,69) GeV and (7 ,gx) = (1.26,0.23), 
the dotted line represents (m x ,m H 0) = (70,60) GeV and ( , y,gx) = (1.22,0.21), and the 
dashed line is (m x ,mf/o) = (60,59) GeV and (7, gx) = (1.23,0.21). The taken values of 
the gauge coupling gx are determined from the observed DM relic density. From the figure, 
we see that when the mass difference m x — m#o becomes larger, due to the boosted H°, 
the flux after the peak of the excess tends to be enhanced and disfavors with the data. 
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Hence, we only focus on m x ~ o. In scheme (b), XX bb is dominant. The result of 
the gamma-ray flux as a function of E x is given in Fig. |6](b) , where gx = 1 is taken, the 
solid and dashed lines stand for (m x ,m H o) = (50,101) and (40,101) GeV, respectively, and 
the value of sin# ~ 0.02 is read from Fig. El(b) for both cases when the observed Dom ^ 2 is 
satisfied. In addition, the value of m H o has been chosen to follow the constraint of the direct 
detection, i.e. m H o ~ 2 m x . For the case of m# o < 2 m x , due to the produced H° being an 
on-shell particle, the annihilation cross section becomes too large to explain the gamma-ray 
excess. Hence, we adopt m H o > 2 m x . 




FIG. 6: Gamma-ray spectrum from dark matter annihilation processes (a) XX —^ H°H° in which 
H° mainly decays into bb and (b) XX —^ bb, where the former corresponds to scheme (a) and 
the latter is scheme (b). The values of slope index 7 are taken as 1.26 [1.22] for (m x ,m^ 0) = 


(70,69) 


Ref. 


(70,60) and (60,59)] GeV and 1.33 for m x = 50(40) GeV. The data are quoted from 


241 ] with ip = 5 degrees . 


Finally, we make some comparisons with the study in Ref. 32] where the stable DM 
candidates are dictated by the custodial symmetry [59fl. Since the trilinear couplings of 


gauge bosons exist in the model given by Ref. 


591 ], besides the annihilation processes which 


we only have in our model, there are also semi-annihilation processes in Refs. ]32|, [59|. With 
the taken values of parameters and the best-fit approach, the authors of Ref. 32j have found 
that the gamma-ray excess is dominated by the semi-annihilation. As a result, DM with 
its mass around 39 — 76 GeV could fit the measured gamma-ray spectrum of the Galactic 
Center. However, the resulted (crw re i) is a factor of 2-3 larger than that of the observed 
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^dm^ 2 - In our approach, with the selected values of m x , e.g. m x = (70,60) GeV in scheme 
(a) and m x = (50,40) GeV in scheme (b), we first constrain the free parameters by using the 
observed Ddm^- and the upper limit of the DM direct detection. With the allowed values of 
parameters, we subsequently estimate the gamma-ray spectrum from the DM annihilation. 
Although the best-fit approach is not adopted in the analysis, our results from the on-shell 
H° production in scheme (a) and m^o ~ 2 m x in scheme (b) are morphologically consistent 
with the gamma-ray spectrum of the Galactic Center. 

In summary, to interpret the excess of the gamma-ray through the DM annihilation, we 
have studied the DM model in the framework of SU(2) X gauge symmetry. To break the 
gauge symmetry, we have used one quadruplet of SU{2) x ■ As a result, the remnant Z 3 
symmetry of SU{2) X leads to the stable DMs, which are the gauge bosons of SU(2) x . Due 
to the mixture of the quadruplet and SM Higgs doublet in the scalar potential, the DM 
annihilation to SM particles is through the Higgs portal. When the observed relic density 
of DM and the limit of the DM direct detection are both satisfied, we find that m x < m\y 
could give a correct pattern for the gamma-ray spectrum. For more specific numerical 
studies, we classify the values of parameters to be scheme (a) with ( m x , m H o) = (70, 69) and 
(60, 59) GeV and scheme (b) with (m x , m H o) = (50,101) and (40, 81) GeV. We show that for 
matching the gamma-ray excess, in scheme (a) it is better to take m x ~ o. If m x — m H o is 
increasing, due to the boosted H°, the gamma-ray flux at the photon energy over the peak of 
the gamma-ray spectrum is enhanced and the resulted flux tends to be away from the data. 
In scheme (b), for avoiding the constraint from the DM direct detection and the production 
of the on-shell H° which causes too large cross section, the condition of m^o > 2 m x is 
adopted. Based on our current analysis, we see that the results of scheme (a) fit the data well. 
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